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Abstract We present a detailed theoretical and ex- 
perimental study on the optical control of a trapped-ion 
qubit subject to thermally induced fluctuations of the 
Rabi frequency. The coupling fluctuations are caused by 
thermal excitation on three harmonic oscillator modes. 
We develop an effective Maxwell-Boltzmann theory which 
leads to a replacement of several quantized oscillator 
modes by an effective continuous probability distribution 
function for the Rabi frequency. The model is experimen- 
tally verified for driving the quadrupole transition with 
resonant square pulses. This allows for the determination 
of the ion temperature with an accuracy of better than 
2% of the temperature pertaining to the Doppler cooling 
limit Tjj over a range from O.STq to bTo- The theory 
is then applied successfully to model experimental data 
for rapid adiabatic passage (RAP) pulses. We apply the 
model and the obtained experimental parameters to elu- 
cidate the robustness and efficiency of the RAP process 
by means of numerical simulations. 



1 Introduction 

Ions trapped in linear Paul traps allow for near-perfect 
realizations of qubits, which can be coherently controlled 
by driving on long-lived transitions between internal elec- 
tronic states with lasers. Important applications range 
from quantum information and simulation experiments, 
quantum metrology or fundamental quantum optics. One 
of the most prominent challenges is given by the require- 
ment to keep one or more ions sufficient well localized 
that sufficiently high fidelities of coherent control opera- 
tions are attained. This is usually achieved by sideband 
cooling of the ions close to the motional ground state of 
one or more vibrational modes. For larger ion numbers, 
however, the number of motional modes increases such 
that complete ground state cooling becomes inconve- 
nient. Furthermore, with the advent of segmented micro- 
traps, it became possible to store an increased number of 



qubits and shuttle the qubits between different trap sites. 
This advantage is bought at the price of heating rates 
which can be larger by several orders of magnitude com- 
pared to conventional linear Paul traps. Furthermore, 
motional energy can be transferred during the shuttling 
operations. Thus, one has to resort to coherent control 
operations which are inherently robust against residual 
motion of the qubits. In order to characterize the depen- 
dence of the fidelity of these operations and find possible 
ways for improvement, a theoretical tool is required for 
simulating the coherent dynamics of the qubit under the 
consideration of multiple excited motional modes. The 
motion leads to fluctuations of the Rabi frequency, which 
represent an inhomogeneous broadening effect when a 
single experimental run is considered, however as ion 
trap experiments are generally conducted by repeating 
a given control sequence and averaging over the qubit 
readout result, the fluctuations manifest themselves as 
an effective homogeneous broadening. We have devel- 
oped an effective Maxwell-Boltzmann theory which en- 
ables us to precisely and efficiently calculate the aver- 
aged dynamics of such a system subject to control laser 
pulses which are shaped in amplitude and frequency. 
This theoretical tool might find further applications for 
closed loop quantum control experiments with shaped 
laser pulses [1] , especially for optimized entangling gates 
working in the thermal regime [2]. 

We present a case study where we use this theoreti- 
cal approach to model electron shelving dynamics on a 
quadrupole transition for a Doppler cooled ion. This has 
direct experimental relevance, as we encode the qubit 
information in the Zeeman split sublevels of the S1/2 
groundstate of a '''"Ca"'" ion [5], and we employ electron 
shelving on the electric quadrupole transition S1/2 -^^5/2 
near 729 nm with a decay time of 1.17 s [3] for qubit 
readout. Population from only one of the qubit levels is 
to be transferred to the metastable state, such that the 
occurrence of resonance fluorescence on a cycling transi- 
tion is conditional on the qubit state. In order to realize 
a transfer process which is sufficiently robust against ex- 
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perimental parameter fluctuations, we make use of the 
rapid adiabatic transfer (RAP) process, where the trans- 
fer pulses are adiabatically switched on and off and the 
laser frequency is swept (chirped) across the atomic res- 
onance [5 . This makes the transfer process inherently 
robust against fluctuations and drifts of the laser fre- 
quency and power, as well as fluctuations of the Rabi 
frequency arising from thermal motion of the trapped 
ion. For the RAP readout, this method offers the possi- 
bility to make a determined trade between robustness, 
efficiency and duration of the process for a given amount 
of maximum available laser power. This manuscript is or- 
ganized as follows: The Maxwell-Boltzmann theory for 
the laser driven qubit is developed in Sec. |2] We then 
give a brief account on our experimental apparatus and 
the used techniques in Sec. [3j In the following Sec. |4j we 
use the Maxwell-Boltzmann theory to accurately model 
Rabi oscillations on the quadrupole transition which are 
subject to thermally induced dephasing. The tempera- 
ture is inferred from the dephasing behavior. We make 
use of the model and the obtained temperature to re- 
produce experimental data for different RAP parameters 
(amplitude and chirp range) without free fit parameters. 
This allows us to use the model for the detailed numeri- 
cal studies presented in Sec. [5j where we investigate the 
dependence of the robustness and maximum fidelity of 
the RAP process on the control parameters. 



2 Effective Maxwell-Boltzmann theory for a 
laser driven qubit 

The interaction picture Hamilton operator of a coher- 
ently and resonantly driven qubit, which is harmonically 
confined in a trap, is given by 

Hi = \hfio&+ (l - ^ {To) a + l)+0 (ry*)) +h.c., (1) 

where J7o is the bare Rabi frequency, and a are the 
raising and lowering operators for the harmonic oscilla- 
tor describing the motion in the trap at frequency lo, 
and rj — cos{a)kyJh/ {2muj) is the Lamb-Dicke factor. 
Here, k is the wavenumber of the driving laser beam 
which makes an angle of a with the oscillation direction 
and m is the ion mass. The Hamiltonian Eq. [l] is a sec- 
ond order expansion in terms of 77, and terms describ- 
ing off-resonant excitation of motional sidebands have 
been neglected, which is justified for rj^/h < 1. Depend- 
ing on the laser beam geometry, the laser can couple to 
three oscillation modes with different frequencies, Lamb- 
Dicke factors and population distributions. This can be 
accounted for in the Hamiltonian Eq. [l] by replacing the 
77^ term by a summation over several modes. If the ion 
is initialized in the ground state \S) and exposed to a 
resonant square pulse of duration t, the probability of 
transferring the ion to the metastable state \D) is given 
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Fig. 1 (color online) EfTective Rabi frequency distribution; 
The symbols show sample probability densities obtained from 
Eq. W\ and the solid lines result from a fit of this data to 
Eq. Is] Data for different temperatures is shown, where the 
temperature is given in terms of the Doppler cooling limit 
Td, see text. The data set for T ~ 2 ■ To (red) corresponds 
to experimental parameters which accurately reproduce the 
measurement results in Fig. [2] The inset shows the linear 
relation between and the ion temperature T. 

by averaging the unitary dynamics over the vibrational 
modes: 

Poit) - (I[Pth{n^,n,) J i(l - cos (^2{„,}t))) , 

(2) 

where the index «=1,2,3 runs over the motional modes, 
such that a threefold summation over the vibrational 
quantum numbers ni=Q..n^^^ has to be carried out. If 
we assume that the vibrational modes are thermally oc- 
cupied, the joint probability to find the ion with a set 
of quantum numbers {rii} is given by the product of the 
three respective Boltzmann distributions 

n" 

Vth{n,n) ^ (3) 

with the mean phonon numbers n" = ksT / {huJi) per- 
taining to the i-th mode oscillating at the frequency LOi. 
The Rabi frequency 

^{n.}=^On*C^(^')' (4) 
i 

is determined by the matrix elements [6] 

^er^^'^ll^i^). (5) 

where the "qi are the Lamb-Dicke parameters for mode 
i and L^{rf') is an associated Laguerre polynomial. The 
triple sum in Eq. [2]is numerically very inconvenient, es- 
pecially if extensive simulations to arbitrary excitation 
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pulse shapes or fitting of experimental data for ion ther- 
mometry are to be performed. Even more, in the case 
of multi-ion crystals, since the number of relevant vi- 
brational modes increases linearly with the ion num- 
ber, the summation approach becomes completely in- 
tractable. The problem of the time evolution on the car- 
rier for a multi-mode thermal distribution has already 
been addressed in Ref. f7|, however, the approach is lim- 
ited to resonant square excitation pulses. Thus, we use 
the alternative approach of averaging over a probability 
distribution function of Rabi frequencies. 
The probability for finding a specific Rabi frequency is 
given by: 



(6) 



Practically, w{^2) is characterized by a set of Rabi fre- 
quencies fi{ni} and their corresponding probabilities. These 
can be evaluated numerically where the motional quan- 
tum numbers are taken into account up to a trunca- 
tion limit. The inaccuracy of this truncation can be es- 
timated by summing the calculated probabilities, which 
should lead to a value slightly smaller than unity. We 
now assume that the differences between the contribut- 
ing Rabi frequencies are smaller than the inverse inter- 
action time and the individual probabilities are small 
due to the many contributing frequencies. Then, the fre- 
quencies can not be individually discerned and a descrip- 
tion by a smooth continuous Rabi frequency probability 
distribution is justified. The probability for attaining a 
Rabi frequency of H during one individual measurement 
is then given by convolution of w{f2) with a Gaussian 
smoothing function: 



1 



V2t 



(7) 



Under the condition that the smoothing parameter a is 
chosen to be larger than the average spacing between 
two Rabi frequencies and smaller than Rabi frequency 
differences that are resolved on experimental timescales, 
'w{i2) is empirically found to be well described by 



V 62„ ) 



(8) 



with the normalization factor A/", which is to be deter- 
mined numerically for given parameters &, i?o- The pa- 
rameter h directly determines the temperature T under 
the assumption that the three motional modes are in 
thermal equilibrium. The relation between T and h is 
determined empirically by calculating for a set of 

fi values and for various temperatures as shown in Fig. 
[TJ We obtain the relation 

where c is uniquely determined by the specific set of 
Lamb-Dicke factors, i.e. by the laser geometry and the 



trap frequencies. The fact that the model function Eq. 
|8] does not perfectly match the actual probability dis- 
tribution, systematic errors in the determination of the 
temperature occur. This leads to uncertainties in the de- 
termination of h and c. From the numerical data shown 
in Fig. [l] and for our specific set of Lamb-Dicke param- 
eters given below in Sec. [s] we find c =4.0(1)-10^. An 
inaccuracy of the temperature of better than 0.02 T^i 
over a temperature range from 0.5 T/j to 5 T^j can be 
claimed. 

The analogy to the Maxwell-Boltzmann velocity distri- 
bution as used in kinetic gas theory becomes clear from 
the distributions shown in Fig. [T] If more than one mo- 
tional degree of freedom is taken into account, a larger 
density of states favors higher thermal energies, whereas 
the ID Boltzmann distribution is monotonically decreas- 
ing with respect to energy, leading to a peaked prob- 
ability distribution. The excess energy translates into 
a Rabi frequency reduction due to an average Doppler 
shift, broadening the optical transition. Quantitative dif- 
ferences to kinetic gas theory arise from the facts that 
i) one is dealing with a harmonically confined instead 
of a free particle, ii) the spatial extent of the thermal 
motion is comparable to the wavelength of the driving 
laser field (in the thermal regime), leading to nonlinear 
coupling matrix elements in Eq. [sj and iii) the confine- 
ment anisotropy and the laser geometry lead to differ- 
ent coupling matrix elements for the vibrational modes. 
The great simplification is that instead of six parame- 
ters r]i,ni, the thermal motion is characterized by only 
one parameter, and the three- fold summation in Eq. [6]is 
replaced by a single integral: 



P^{t)^- [ " Wi,{f2){l-cos{Qt))dn. (10) 
2 Jo 



The relation between the bare Rabi frequency i7o and 
the experimentally determined pulse time when the first 
excitation maximum occurs Tmax is given by: 



(11) 



as obtained from dwb/dn\^ = 0, corresponding to the 
effective Rabi frequency D — 7r/rmax with the largest 
probability. The practical application of our model is as 
follows: 

For a given laser-ion interaction setting (beam geome- 
try, ion number, trap frequencies), one has to calculate 
the data characterizing the w{n) Eq.[6]for a set of exam- 
ple temperatures as shown in Fig. [l] Then, the smoothed 
probability distribution Eq.[8]is to be determined at sam- 
ple points such that the b parameter can be extracted 
from a nonlinear regression. Once the relation between 
T and b is established and the bare Rabi frequency is 



calculated from Eq. 11 the model can be applied e.g. as 



in Eq. 10 to simple Rabi oscillations. 
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3 Experimental Apparatus 

In our experiment, we Iceep a single ^*'Ca+ ion in a 
microstructured linear Paul trap l8:,,_9j. The oscillation 
frequency pertaining to the axial mode of vibration is 
Wax = 27r-1.35 MHz, while the radial modes oscillate at 
Wiad = 27r • {2.4,3.0} MHz. A quantizing magnetic field 
at 45° to the trap axis provides a Zeeman splitting of 
Az «27r-18 MHz between the mj = ±1/2 levels of the 
Si/2 ground state, which are henceforth denoted as ||) 
and It)- The ion is Doppler cooled on the S1/2 — ^Pi/2 
cycling transition near 397 nm, while resonance fluores- 
cence from this transition is monitored both on an EM- 
CCD camera and a photomultiplier tube (PMT). Co- 
herent dynamics are driven on the 5'i/2 ^5/2 E2 
quadrupole transition near 729 nm, where laser light 
from an amplified diode laser system is switched and 
modulated by a double-pass acousto-optical modulator 
(AOM) running at frequencies around 80 MHz. Up to 
100 mW of laser light is focused onto the trapping site 
with a lens of 200 mm focal length, resulting in spot size 
at the focus of about 30 ^m FWHM. The propagation 
axis is aligned to be at 45° to the trap axis and orthogo- 
nal to the quantizing magnetic field. This leads to Lamb- 
Dicke factors for the three vibrational modes of \r\i\ = 
{0.059,0.031,0.028}. The polarization of this laser is cho- 
sen to be at an angle of 45° to the plane defined by the 
propagation direction and the magnetic field. According 
to the selection rules for E2 transitions, this allows for 
driving of the transitions |t) — )■ I-D5/2, mj = +5/2) and 
It) \Dz/2,mj = +3/2). The AOM is supplied with 
r.f. signals from a DDS-based synthesizer which en- 
ables shaping of the laser pulses in frequency, phase and 
amplitude. 

The experiments are carried out as follows: the ion is 
Doppler cooled on the cycling transition and initialized 
in It) by optical pumping on the same transition. Then, 
the 729 nm excitation pulse is irradiated onto the ion, af- 
ter which the readout is performed by irradiation on the 
cycling transition. If the ion was excited to the metastable 
state, fluorescence rates of typically about 30-10'^ counts/s 
are observed, whereas background count rates of about 
4-10^ counts/s are measured if the ion resides in the 
metastable state. After the readout, the ion is reset by 
quenching the metastable state by irradiation of laser 
light near 854 nm driving the D5/2 -^^3/2 dipole tran- 
sition. For a given set of excitation pulse parameters, 
the sequence is repeated 200 times such that the excita- 
tion probability is determined with a shot noise limited 
inaccuracy of at maximum 3.5%. 



4 Experimental results 

The measured excitation probability on the transition 
It) — >■ |-D5/2,TOj = +5/2) versus pulse duration for square 
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Fig. 2 (color online) Coherent dynamics on the quadrupole 
transition: The fraction of population in the metastable 
state is plotted against the duration of a square excita- 
tion pulse with the laser frequency tuned to the |t) 
ID5/2, wij = +5/2) transition. One can clearly see a very 
rapid dephasing, which is due to the interaction with three 
thermally excited vibrational modes. Note that the dephas- 
ing behavior, i.e. the envelope of the oscillations, cannot be 
described by a Gaussian or exponential decay. 

pulses is shown in Fig. [2] along with the fit to the model 
Eq. [To] The period allows for a direct extraction of r„iax 
= 4.93(5) /is and therefore allows for the determination 



of the bare Rabi frequency according to Eq. 11 which is 
found to be fl =27r-105(l) kHz for the shown data. With 
this knowledge, the parameter 6, characterizing the dis- 
tribution of the Rabi frequencies, can be fixed from a fit 
of the signal to the model Eq. (To] We find b =7.1-10"'', 
which is consistent with T « 2.0 • To if we assume that 
all modes are cooled to the same temperature, leading 
to mean phonon numbers of {rii} = {16.9,9.5,7.6} . The 
temperature is found to be two times larger than the 
theoretical Doppler cooling limit of Td wO.55 mK for an 
ideal two-level system. This discrepancy is attributed to 
a competing micromotion-induced Doppler heating ef- 
fect. Further reasons for decreased Doppler cooling efh- 
ciency are the invalidity of the two-level approximation 
due to the presence of the Zeeman splitting (which is 
comparable to the linewidth of the cooling transition) 
and the additional decay channel to the D3/2 state. 
For the detailed evaluation of the RAP pulses, we first 
gauged the bare Rabi frequency f2o versus the precisely 
controllable r.f. voltage which is supplied to the AOM. 
The maximum excitation time is measured by recording 
the first oscillation period of a signal such as the one 
shown in Fig. [2] for a set of different r.f. amplitudes. J7o 
is then determined by extracting the effective Rabi fre- 
quency as explained above, and making use of Eq. 
[TT| The resulting Rabi frequencies are used for a fit to 
a third order polynomial function, such that the Rabi 
frequency for any r.f. amplitude can be obtained. In the 
experiment, we use RAP pulses with temporally varying 



^ VFG-150, Toptica Photonics, Graefelfing, Germany 



Interaction of a Laser with a Qubit in Thermal Motion and its AppUcation to Robust and Efficient Readout 



5 




100 200 300 400 500 



Peak Rabi frequency (angular) [kHz] 

Fig. 3 (color online) RAP efficiency versus peak Rabi fre- 
quency for different chirp ranges: The plot shows resulting 
transfer efficiencies while the peak amplitude is scanned. 
The solid lines are obtained from a numerical solution of 
the time-dependent Schodinger equation for the qubit, in- 
cluding thermal effects by averaging over a distribution Rabi 
frequencies (see text). For small chirp ranges, one observes a 
Rabi oscillation-like behavior, while adiabaticity is attained 
for chirp ranges of fOO kHz or larger. For even larger chirp 
ranges, the increase in robustness is bought at the expense 
of a higher power requirement. Note that no free parameters 
were used for the simulation, all parameters were inferred 
from the pulse width scan measurement of Fig. [2] and the 
power gauge measurement. 
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Fig. 4 (color online) Parasitic shelving: The figure shows 
the shelved population versus peak Rabi frequency for the 
same parameters as in Fig. |3] The central frequency is res- 
onant with the It) — >■ \D^/2,mj — -1-5/2) transition, but the 
ion was initialized in |^) such that no population should be 
transferred under ideal conditions. Data for different chirp 
ranges is shown. It can be seen that population transfer is 
insensitive to the chirp range, which suggests a completely 
incoherent transfer mechanism. 

eraging over the thermal Rabi frequency distribution Eq. 
[sjis done. The parameter b =7.1T0~'* is chosen accord- 
ing to measurement results from Fig. [2] such that no free 
fit parameter is occurring. 



Rabi frequency Q{t) and detuning 5{t) according to 



n{t) 

5{t) = nrct/T„, 



(12) 



where q'^'^^^ the peak amplitude known from the power 
calibration, determines the duration of the pulse, which 
has a Gaussian envelope, and the linear frequency chirp 
across the resonance is determined by the chirp range r^. 
The pulses are truncated in time at i±2Tcr, and both fre- 
quency and amplitude are changed in 50 discrete time 
steps. Fig. |3] shows the excitation probability, i.e. the 
probability to find the ion in the metastable state, ver- 
sus peak Rabi frequency for different chirp ranges and 
fixed T„ =50 ^s, resulting in a 4 /xs sample duration. It 
is important to note that a phase-continuous frequency 
switching mode of the r.f. synthesizer is used, where the 
phase of each r.f. sample is chosen such that the output 
signal is continuous. One clearly observes the transition 
from a Rabi-oscillation behavior to a monotonic increase 
in excitation with the amplitude as the chirp range be- 
comes larger and adiabatic following conditions are at- 
tained at Tc =100 kHz. For larger chirp ranges, higher 
amplitudes are needed to attain maximum excitation, as 
the effective time on resonance is decreased. Fig. [3] also 
displays the results of simulations for each chirp range, 
where the time-dependent Schrodinger equation is solved 



numerically for the excitation pulse Eq. 12 while the av- 



The fidelity of readout process can also be deterio- 
rated if population is excited from the wrong qubit state, 
i.e. 14,) in our case. Conceivable mechanisms for this 
transfer are off-resonant excitation on parasitic transi- 
tions, resonant excitation of motional sidebands of these 
transitions, or resonant excitations by frequency com- 
ponents of the excitations pulse which arise from the 
discrete sampling. Fig. |4] shows the population trans- 
ferred to the metastable D5/2 state by similar pulses as 
in Fig. [sj but for initialization in the jj,) level. One can 
see that an error of about 1% is present at a peak Rabi 
frequency of 200 kHz, and the excitation is independent 
of the chirp range. The transition connecting to jj,) with 
the lowest frequency offset is the ||) — > I-D5/2, "ij + 3/2) 
transition, which is separated from the driven jf) — >■ 
|£'5/2,m-j = +5/2) by about 8 MHz. This is far beyond 
the used chirp ranges, i.e. out of the bandwidth of the 
excitation pulses. The simulation yields only population 
transfer values up to 10^^ for the maximum amplitude, 
where the values are strongly increasing with the chirp 
range. Thus, the excitation mechanisms via resonant fre- 
quency components and resonant excitation of motional 
sidebands can therefore be excluded. It can therefore 
be concluded that the parasitic transition is excited by 
incoherent amplified spontaneous emission (ASE) back- 
ground from the amplified laser. Thus, this is a purely 
technical artifact which can be overcome by employing 
a filter cavity. 
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Fig. 5 (color online) Simulation of the RAP pulse effect for systematic deviation of experimental parameters. The plots 
show simulation data of the logarithmic infidelity logj^Q (1 — poite)), i.e. a measure for the population not transferred to the 
metastable state. The infidelities are plotted versus a factor describing static errors of the Rabi frequency, y — Qq/H'i^'^^^ 
with respect to the one obtained from the calibration q}^^^'^ , and versus an extra detuning 5' describing an error of the laser 
frequency with respect to the atomic resonance, which is given in units of the chirp range of 27r-100 kHz from plot b). Plot 
a) shows the infidelity data for zero chirp range, i.e. in the Rabi oscillation regime. The peak Rabi frequency is chosen to 
be n^"''^'' « 27r-332 kHz. One can clearly see that the parameter regions for which a considerably high population transfer is 
accomplished are rather low, and a minimum infidelity of about -2 is attained. Plot b) shows the infidelity for a chirp range 
of 27r-100 kHz and a peak Rabi frequency of fig'^^^^ « 27r-221 kHz, corresponding to the situation where adiabatic following 
conditions are fulfilled and the transfer efficiency saturates with respect to laser power as inferred from the experimental data 
Fig. Is] Plot c) shows the same simulation for a chirp range of 27r- 150 kHz and a correspondingly increased Rabi frequency of 
n^'^''^^ 27r-332 kHz. 



5 Numerical study of fidelity and robustness 



pulse is averaged according to 



The RAP pulses are characterized by the parameters 
duration, amplitude and chirp range, and a quantitative 
understanding of the role of these parameters is required 
to achieve predetermined fidelity and robustness val- 
ues. As the situation is complicated due to the presence 
of thermal fluctuations, the Maxwell-Boltzmann model 
proves to be a useful tool for a precise numerical char- 
acterization of the dynamics. For a given thermal Rabi 
frequency reduction factor x, we numerically solve the 
Schrodinger equation for the qubit in ODE form: 



(13) 



(14) 



c(-) = i((^(t)+y)cf +xyr2(i)c(-) 
(Sit) + S')4r^ + xyf2{t)c'^J'^ 



Here, Cg{t) and Ce{t) are the amplitudes of the ground 
and excited state, respectively, and the system is as- 
sumed to be perfectly prepared in the electronic ground 
state prior to the excitation pulse. S' represents an addi- 
tional offset detuning arising from possible drifts of the 
laser frequency, and y is an additional constant scaling 
factor for the Rabi frequency from laser intensity drifts. 



For a given parameter set S',y, Eqs. 13 are numerically 
solved for < a; < 1 in steps of Ax 
suiting transfer efficiency p^^\te) = 



= 0.01 and the re- 
i''\te)\^ after the 



with Wb{f^) from Eq. [s] 

Fig. [5] shows the resulting logarithmic infidelity, i.e. the 
population residing in the initial state after the transfer 
pulse, versus the frequency and amplitude error param- 
eters for three different chirp ranges. The duration of 
T(j =50 /iS, the b parameter and the sample number are 
chosen to be constant and identical to the experimental 
parameters from Fig. [3j Dephasing Rabi oscillation be- 
havior is seen for zero chirp range, and it becomes clear 
that the frequency and amplitude have to be rather pre- 
cisely calibrated and stabilized in order to keep the trans- 
fer efficiency sufficiently high. By contrast, chirp ranges 
of 100 kHz and more lead to a very robust behavior. For 
an increased chirp range of 150 kHz, extra laser power 
is needed to keep adiabatic following conditions satis- 
fied, but the maximum transfer efficiency is increased 
by more than one order of magnitude, and the robust- 
ness window is greatly enhanced. This illustrates how the 
RAP process allows to make use of extra laser power to 
obtain an increased robustness and a higher maximum 
efficiency. Note that transfer infidelities below 0.1% can 
hardly be actually observed in the experiment as other 
limitations from the state preparation and fluorescence 
readout steps obscure the result. 
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6 Conclusion and Outlook 



We have developed a Maxwell-Boltzmann theory for ther- 
mally induced Rabi frequency fluctuations in the coher- 
ent excitation of a single trapped ion with a laser. The 
model is successfully applied to reproduce the dephas- 
ing behavior of Rabi oscillations and to RAP excitation 
pulses with variable parameters. The latter opens the 
possibility to make the right choice of pulse parameters 
to achieve well defined behavior in terms of robustness 
and efficiency if the ion temperature and the relation be- 
tween Rabi frequency and laser power is known. Based 
on the model, we have carried out numerical simula- 
tions which allow for the determination of the robustness 
against experimental parameter fluctuations. This makes 
the RAP pulse a very reliable tool for population trans- 
fer in quantum experiments with trapped ions. Further- 
more, the model offers a new thermometry scheme which 
operates from closely above the Lamb-Dicke regime to 
temperatures several times larger than the Doppler cool- 
ing limit temperature. It closes a gap between two estab- 
lished methods: The fully coherent measurement method 
which compares the Rabi frequencies of the red and blue 
motional sidebands on a given motional mode, which 
works only in the Lamb-Dicke regime and requires that 
either the interaction with spectator modes is suppressed 
or that these are also cooled deeply into the Lamb-Dicke 
regime. A simpler scheme based on time-resolvcd fluo- 
rescence readout [TU1[TT| does not require the ability of 
coherently driving long-lived transitions, but ceases to 
yield precise answers at temperatures close to and be- 
low the Doppler limit. 

Our technique also opens the possibility for the appli- 
cation of quantum control strategies such as optimal 
control theory for thermal ensembles, e.g. for improv- 
ing the robustness of quantum gates. Entangling gates 
are usually the operations in quantum experiments with 
trapped ions which are the most difficult to realize with 
considerable fidelity, as the interplay between internal 
and external degrees of freedom requires the precise knowl- 
edge and control of many parameters. Residual thermal 
excitation of the motional modes leads to a loss of con- 
trollability of the system • A recent example of how 
coherent control strategies can cure the effect of control 
parameter fluctuations for an entangling quantum gate 
is given by |13| . The averaging of the coherent dynam- 
ics as in Eq. [lOjfor general Hamiltonians in conjunction 
with optimal control theory [M] could lead to the devel- 
opment of gate schemes which are more robust in the 
thermal regime than current ones 
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